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By utilizing the property of the supersymmetric structure in the two-level multiphoton Jaynes- 
Cummings model, an invariant is constructed in terms of the supersymmetric generators by work- 
ing in the sub-Hilbert-space corresponding to a particular eigenvalue of the conserved super- 
' symmetric generators. We obtain the exact solutions of the time-dependent Schrodinger equa- 

I tion which describes the time-dependent supersymmetric two-level three-photon Jaynes-Cummings 

. model (TLTJCM) by using the invariant-related unitary transformation formulation. The case under 

the adiabatic approximation is also discussed. 
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I. INTRODUCTION 

^ ^ The interaction between a two-level atom and a quantized single-mode electromagnetic field can be described by the 
1 Jaynes-Cummings model (JCM) [1] which has been applied to investigate many quantum effects such as the quantum 
' collapses and revivals of the atomic inversion, photon antibunching, squeezing of the radiation field, inversionless 
I light amplification, electromagnetic induced transparency [2-5], etc.. In addition to the standard JCM, there exists 
' another type of JCM which possesses supersymmetric structure [6,7]. In this generalization of the J-C model, the 
. atomic transitions are mediated by k photons. Singh has shown that this model can be used to study multiple 
^ atom scattering of radiation and multiphoton emission, absorption, and laser processes [8]. Some authors introduced 
I a supersymmetric unitary transformation to diagonalize the Hamiltonian of this supersymmetric JCM and obtain 
] the eigenfunctions of the stationary Schrodinger equation [9, 10]. It is of great interest to investigate the geometric 
I . phase factor of the time-dependent JCM since the geometric phase factor [11] appears only in systems with the time- 
' dependent Hamiltonian. It can be verified as shown in the present paper that the exact solutions and the geometric 
. phase factor of the two-level JCM whose Hamiltonian has time-dependent parameters can also be obtained by making 
' use of the generalized invariant theory [12]. For simplicity and convenience, in this paper we only investigate the 
• • I time-dependent supersymmetric two-level three-photon Jaynes-Cummings model (TLTJCM). This method can also 
. ^ ' be generalized to the time-dependent JCM with more than three photons. 

. The invariant theory which is appropriate for treating time-dependent systems was first proposed by Lewis and 
?H ' Riesenfeld (L-R) [13] in 1969. In 1991, Gao et al. generalized the L-R invariant theory and proposed the invariant- 
related unitary transformation formulation [12,14]. In this formulation the eigenstates of the time-dependent invariants 
are replaced with those of the time-independent invariants through the unitary transformation and the exact solutions 
(which contain the dynamical and geometric phase factor) of the time-dependent Schrodinger equation are obtained. 
Many works have shown that the invariant-related unitary transformation approach is a powerful tool for treating 
time-dependent problems and geometric phase factor [15-17]. 

II. THE EXACT SOLUTIONS OF THE TIME-DEPENDENT TLTJCM 

The Hamiltonian of the TLTJCM under the rotating wave approximation is given by 

H{t) = u;it)a^a + + g{t){a^fa^ + 5* (2.1) 



*E-mail address: jqshen®coer. zju.edu.cn 



1 



where a'^ and a are the creation and annihilation operators for the electromagnetic field, and obey the commutation 
relation [a, a^] = 1; a± and CTz denote the two-level atom operators which satisfy the commutation relation [iJ2, (j±] = 
±2f7± ; g(t) and g*{t) are the coupling coefficients and 3 is the photon number in each atom transition process; 
ujoit) and uj{t) are respectively the transition frequency and the mode frequency. All the parameters in Eq. (2.1) are 
time-dependent, and the time-dependent Schrodinger equation for this system is 

i^-^ = m\m)s- (2-2) 

The supersymmetric structure can be found in the TLTJCM by defining the following supersymmetric transforma- 
tion generators [9,10]: 

«=(«')'"--( (a?). 1!)' «'=''V=(2f). P.3) 

It is easily verified that {N, N , Q, Q^) form supersymmetric generators and have supersymmetric Lie algebra proper- 
ties, i.e., 

g2 = (Qt)2 ^ 0, [Q\ Q] =N'a,, [n, n'] = 0, [N, Q] = Q, 
[7V,Qt]=_Qt, {q\q}=n', {Q,a,} = {Q\a,}=0, 

[Q,a,]=2Q, [Qla,]=-2Q\ (gt _ Q)^ = _7v', (2.4) 

where {} denotes the anticommuting bracket. By the aid of Eqs. (2.3) and (2.4), the Hamiltonian (2.1) of the 
TLTJCM can be rewritten as 

Hit) = u,m + Mt)-Mt) ^^ ^ g^^^Q ^ g.^^^Q, _ a|) _ (2.5) 

According to the L-R invariant theory [13], one should first construct an invariant I{t) in order to show the solvability 
of Eq. (2.2). A Hermitian operator I{t) is called invariant if it satisfies the following invariant equation 

^ + l[I{t),H{t)]=Q, (2.6) 

and the eigenvalue equation of the time-dependent invariant is given by 

/(i)|A„,t) = A„|A„,t), (2.7) 

where ^ = 0. It can be seen from the invariant equation (2.6) that I{t) is the linear combination of N,az,Q and 
Q^. However, it should be emphasized that the generalized invariant theory [12] can only be applied to study the 
system with the quasialgcbra defined in [18]. Unfortunately, there is no such quasialgcbra for the TLTJCM in Eq. 
(2.4). Many problems have been solved in Quantum Mechanics by working in the sub-Hilbert-space corresponding to 
a particular eigenvalue of the Hamiltonian [19] . We show that in the case of the TLTJCM, a generalized quasialgebra, 
which enables one to obtain the complete set of the exact solutions for the TLTJCM, can also be found by working 
in a sub-Hilbert-space corresponding to a particular eigenvalue of the supersymmetric generator N . 
Using a^{a^Y \m) = ^"'^i^^' |m) and (a'^Ya'^ |to -|- 3) = |m -I- 3) ,one can arrive at 

with Am — ^"^j^P' ■ One thus obtains the supersymmetric quasialgebra {N,Q,Q^ ,az) in the sub-Hilbert-space corre- 
sponding to the particular eigenvalue A^ of iV , by replacing the generator N with A™ in the commutation relations 
of Eq. (2.4), namely, 

[O^Q]=A„a„ {Q^Q} = A„, (Q^-Q)' = -A„. (2.9) 
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In accordance with the invariant theory, we chose the invariant I{t) to be of the form 

I{t) = —[ex.p{-i(l))Q + exp( i<l))Q^]+ cos ea„ (2.10) 

where 9 and (j) are time-dependent parameters. Substitution of the expressions (2.5) and (2.10) for I{t) and H{t) into 
Eq. (2.6) leads to the following set of auxiliary equations 

1 

6cos6exp{—i(j)) — i(psm6 e^p{—i(j)) + — wq) sm6exp{—i(j)) — 2gXm cos6] = 0, 

9 - ixi[ge^p{i(l>) - g* exp(-i(?i)] = 0, (2.11) 

where the dot denotes the time derivative. The two time-parameters 9 and (f) in I{t) are determined by these two 

auxiliary equations. 

Using the invariant-related unitary transformation method [12], we define the unitary transformation operator as 
follows 

V{t) = exp[l3(t)Q- 0*{t)Q^] (2.12) 

with l3*{t) being the complex conjugation of f3{t). With the help of the commutation relations (2.4) and by the 
complicated and lengthy computations, it can be found that if P{t) and f3*{t) satisfy the following equations 

^^_fexpH^^ ^^^JeMict>)^ 

\ 2 \ 2 

then the following time-independent invariant ly can be obtained 

Iv = V\t)I{t)V{t) = a,. (2.14) 
Correspondingly, the Hamiltonian (2.5) can be transformed into 

Hv{t)^V\t)H{t)V{t)-V\t)^^V{t) 

OJ 1 1 

= ioN + —{az - 1) + {--Xm[ge:x.pii(f>) + g* exp{-i(l))] sm9 + 

+ ^^^^^cos^-^(l-cos0)K. (2.15) 
by using the Baker-Campbell-Hausdorff formula [20] 

= + i[|L, L] + 1[[|l, L],LH i[[[|L, L],L],LH- (2.16) 

with V{t) = cxp[iy(t)]. Under this unitary transformation2.12 the time-dependent Schrodinger equation (2.2) is then 
transformed into the following form 

i^^^^=Hyit)M,, (2.17) 

where 

|M/(t)), = y(t)|A„,i),o. (2.18) 

One may show that the particular solution |A„,f)^ of Eq.(2.2) differs from the cigcnfimction |A„,t) of the invariant 
I{t) only by a phase factor exp[z^^(i)]. Then tlic general solution of the Schrodinger equation (2.2) can be written as 

!*(*)). = E eMiMt)] \^n,t), (2.19) 

n 

where 
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•n{t)= I (^Xn,t\^-^-Hit)\x„,t)dt, 



Cn = {Xn,t = 0\^{t = 0))^. (2.20) 

It is easy to verify that the particular solution \Xn,t)^Q of the time-dependent Schrodinger equation (2.17) is different 
from the eigenfunction |A„) of ly only by the same phase factor exp[i(/>„(t)] as that in Eq.(2.19), i.e., 

\Xn,t),Q = exp[i(^„(t)] |A„) . (2.21) 

Substitution of |A„,t)^o of Eq. (2.17) into Eq. (2.21) yields 

-^{t)\Xn)=Hv{t)\Xn), (2.22) 

which means that Hv{t) differs from Iv{t) only by a time-dependent multiplying c-number factor. Then the particular 
solution of Eq. (2.17) can be easily obtained by calculating the phase from Eq. (2.22). 

The eigenstates of corresponding to the eigenvalue a = +1 and a = —1 are (J) and (°), and the eigenstate of 
N' is (ij,"^3)) corresponding to Eq. (2.8). From Eq. (2.20), (2.21), (2.22), we obtain two particular solutions of the 
time-dependent Schrodinger equation of the TLTJCM which can be written in the following forms 



where 



and 



and 



where 



*m,a=+i(t)), = exp{i ^ (i')+ {t)\dt}V{t) {^"^^^ (2.23) 



^d,a=+i {t) = {m+^)w{t) - ^Xi{g{t)exp[i(l>{t)]+g*{t)exp[-i^{t)]}sme{t) 

^ Mt)-Mt) ^^^^^,^ (2.24) 



Vg,a=+i it) = -^[1 - cos0(i')]; (2.25) 



|*m,a=-l(t)), = exp{^ ^ [^4,.=-! {t)+ fg,a=-l {t)]dt}V{t) (^^^ ^ (2.26) 



3 1 ^ 

<a=-i {t) = (m+ -)w(t') + -Xl{g{t)exp[i<l>{t)]+g*{t)ei^p[-i<f>{t)]} sin 9{t) 

_Mt)-Mt)^^^,^^'^ (2.27) 



and 



(t) = ^[1 - cos9it)]. (2.28) 

Generally speaking, in Quantum Mechanics, solution with chronological-product operator (time-order operator) P, 
namely, where the time-evolution operator U{t) = Pexp[i Jq H{t )dt ],is often called the formal solution. In the 
present paper, however, the solution of the Schrodinger equation governing a time-dependent system is sometimes 
termed the explicit solution, for reasons of the fact that it does not involve time-order operator. But, on the other 
hand, by using Lewis-Riesenfeld invariant theory, there always exist time-dependent parameters, for instance, 6 and 
(p in this paper which are determined by the auxiliary equations (2.11). Traditionally, when employed in experimental 
analysis and compared with experimental results, these nonlinear auxiliary equations should be solved often by means 
of numerical computation. In view of the above reasons, the concept of explicit solution is understood in a somewhat 
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relative sense, namely, it can be considered explicit solution when compared with the time-evolution operator involving 
time-order operator; whereas, it cannot be considered completely explicit solution for it is expressed in terms of some 
time-dependent parameters which should be obtained via the auxiliary equations. Hence, conservatively speaking, 
we regard the solution of the time-dependent system presented in the paper as exact solution rather than explicit 
solution. 

These above two particular solutions of the Schrodinger equation (2.2) contain the corresponding dynamical phase 

factor exp[i Vd.a {t )dt ] and geometric phase factor exp[i '^g.a {t )]dt with a = ±1. Apparently, it can be 
seen that the former phase factor depends on the transition frequency Wo(i) and the mode frequency Lu(t) as well as 
the coupling coefficients g{t) and g*{t), whereas the latter is immediately independent of these frequencies and the 
coupling coefficients. 

One of the theoretical applications of the exact solution of time-dependent Schrodinger equation is constructing the 
time-dependent coherent state [21]. For instance, one result in this paper may be given as follows 

I (t)) =exp( 

= exp( 

with 5 being a time-independent parameter and </'(t=-i-i (t) =fd,(7=+i (t)+ Vg,a=+i {t)- Time-dependent coherent 
state is believed to be useful in investigating the classical properties of supersymmetric Jaynes-Cummings model in 
quantum optics. 
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III. DISCUSSION UNDER THE ADIABATIC APPROXIMATION 

The geometric phase factor in the adiabatic evolution of quantum systems was first discussed by Berry [22] . Simon 
then showed that this phase is in connection with a holonomy of the connection in the Hcrmitian line bundle over 
the parameter space [23]. This adiabatic phase, referred to as the Berry phase, has attracted many attentions of 
the investigators in various branches of physics [24]. Here we investigate the cases of the adiabatic limit. Under the 
adiabatic limit, we assume that the time derivative of 6 vanishes, namely, 

^ = 0, (3.1) 

then the following equations can be derived from the auxiliary equations (2.11) 

3 = Iff] exp(-i0), g* = \g\exp{i(j)), {3u; - ujq - (p) sin 9 = 2 \g\ cos 9 . (3.2) 
Inserting Eqs. (3.2) into Eq. (2.10), we obtain 

= [g{t)Q + g*{t)Q^ + _l(3u; - a;o - <^)a.]. (3.3) 

M — UJQ — (j) ^ 

In accordance with the definition of an invariant, i.e., Eq. (2.6), one can draw a conclusion that an invariant can 
be regarded as the Hamiltonian of the adiabatic-evolution system. Making a comparison between Eq. (2.5) and Eq. 
(3.3), one can see that if 

= 0;, (3.4) 
the following relation between H{t) and I(t) can be obtained 

Further analysis shows that for the general three-generator Lie-algebraic systems, the invariant /(i) is just considered 
the Hamiltonian i?(t) which can be expressed as 

H[t) ~ m (3.6) 
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in the adiabatic-evolution process. According to Eq. (2.6) and Eq. (2.7), the generator iV is a time-independent 
invariant, while I{t) is a time-dependent invariant. Since the eigenstate of N can be rewritten as 

'"^^ ^ \m)(!)+\m + 3)(') (3.7) 



m + 3)J \0J \1 

which is apparently not the eigenstate of /(<), even in the adiabatic evolution, the problem of eigenvalue is not very 
necessary as in the stationary Schrodinger equation. An invariant which satisfies Eq. (2.6) is just a conserved operator. 
The product of two invariants, e.g., I{t)N is also an invariant [12]. More invariants can be constructed in terms of 

0{t) = U\t)OU{t) (3.8) 

where O is an ordinary operator and U{t) is the time-evolution operator which is given 

U{t) = V{t) exp{i f[^d,a (*')+ Vg^a {t )]dt } . (3.9) 
^ Jo 

Here the invariant-related unitary transformation formulation provides the evolution operator with an explicit expres- 
sion rather than a formal solution of equation 

i^=H{t)U. (3.10) 

It should be noted that, using Eq. (3.4) and (3.1), the geometric phases Jg '^g^a {t)\dt can be rewritten as 

a /■* 

^g,ait) = --{l-cose) ujdt (3.11) 
^ Jo 

with C7 = ±1. The geometric phases (or Berry phase) in a cycle (i.e., one round trip) is 

(^g,,(T) = -|27r(l-cose), (3.12) 

where 27r(l — cos^^) is the solid angle over the parameter space of the invariant I{t), which unfolds the geometric 
meanings of the phase factor. The right-hand side of Eq. (3.12) is analogous to the magnetic flux produced by a 
monopole of strength — ^ existing at the origin of the parameter space. This, therefore, implies that geometric phase 
differs from dynamical phase and it involves the global and topological properties of the time evolution of a quantum 
system. 



IV. CONCLUDING REMARKS 



In the present paper we have constructed an invariant in the sub-Hilbert-space corresponding to a particular 
eigenvalue of the conserved operator (time-independent invariant) N and obtained the exact solutions of the time- 
dependent supersymmetric TLTJCM by making use of the invariant-related unitary transformation formulation. This 
formulation replaces eigenstates of the time-dependent invariants by those of the time- independent invariants through 
the unitary transformation. In view of the above calculation, we can see that this unitary transformation formulation 
has some useful applications. 

The exact solutions as well as their geometric phase factors of the time-dependent single- and two-photon cases can 
be obtained by using the present method. Since the three-level two-mode Jaynes-Cummings model plays an important 
role in Quantum Optics [25] , the supersymmetric structure and the exact solutions of the time-dependent three- level 
two-mode multiphoton JCM deserves further investigations by the formalism suggested in the present paper. It is also 
interesting and necessary to obtain the exact solutions of the supersymmetric TLTJCM without the rotating wave 
approximation, by using this invariant-related unitary transformation formulation. 
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